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Abstract. In this article we give an account of a method of smoothing spatial inhomo- 
geneous data sets by using wavelet reconstruction on a regular grid in an auxilliary space 
onto which the original data is mapped. In a previous paper by the present authors, we 
devised a method for inferring the velocity potential from the radial component of the 
cosmic velocity field assuming an ideal sampling. Unfortunately the sparseness of the real 
data as well as errors of measurement require us to first smooth the velocity field as ob- 
served on a 3-dimensional support (i.e. the galaxy positions) inhomogeneously distributed 
throughout the sampled volume. The wavelet formalism permits us to introduce a min- 
imal smoothing procedure that is characterized by the variation in size of the smothing 
window function. Moreover the output smoothed radial velocity field can be shown to 
correspond to a well defined theoretical quantity as long as the spatial sampling support 
satisfies certain criteria. We argue also that one should be very cautious when comparing 
the velocity potential derived from such a smoothed radial component of the velocity field 
with related quantities derived from other studies (e.g : of the density field). 
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1. Introduction 

An 'inverse' problem posed frequently in cosmology is to construct the velocity po- 
tential (assuming an irrotational velocity field) from observations of only the 'radial' 
component of the velocity field. We have argued elsewhere [7] that wavelet transforms 
provide a convenient formalism for this procedure, but subsequently in the course of our 
practical investigations we found that the signal is easily distorted by that due to the 
inhomogeneity of a non-uniformly sampled data set. Such a sample is constituted by a 
set of galaxy peculiar velocities distributed (as they in practice always are) over a nearly 
random support, and it is to this example that we address most of our remarks in this 
paper. However the solution that we propose, namely a non-linear algorithmic mapping 
from the real space into a fictitious but uniformly sampled space may be of general interest 
and is therefore reported here. 

The ultimate cosmo logical objective is to extract information concerning the fluctu- 
ations of the total mass in the universe <5t(x) from the observed radial component iv(x) 
of the cosmic peculiar velocity field v(x), assuming that the latter has been created by 
fluctuations in the self-gravitating matter. One of the goals is indeed, by comparing the 
fluctuations of the total mass obtained as above with the fluctuations of the luminous 
mass <5;(x) obtained from the study of the distribution of the galaxies in the universe, to 
characterize the distribution of the dark matter in the universe <5^(x) = 5*(x) — ^(x) . 

According to the standard models of large scale structure formation, there are good 
reasons to believe that the cosmic velocity field is irrotational at least above some ill- 
defined scale s c . Bertschinger and Dekel [1] have pointed out that consequently above the 
scale s c , the cosmic velocity field v(x) may be derived from a potential (v(x) = V ^(x), 
i.e : v is curl-free), and that this kinematical potential $ can be extracted by integrating 
the observed radial component of the velocity field tv(x) along the line-of-sight : 



Moreover, above the scale s c , the velocity fluctuations still belong to the linear regime and 
the total mass density field S t (x) is thus linked to the kinematical potential $(x) through 
the Poisson equation (<5 t (x) oc V 2 $(x)). 

In our previous paper [7] , we devised a method based on the properties of the wavelet 
transform for inferring the kinematical potential from the observed radial component of 
the cosmic velocity field. We summarize in section 2 why our method is well adapted for 




(1.1) 
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solving the problem of the choice of the scale s c ,and how we can directly derive from the 
radial velocity field quantities such as the laplacian of the kinematical potential (and thus 
the fluctuation in the total mass). 

Unfortunately, as was observed above, the radial velocity field is measured with large 
statistical uncertainties and is sampled on a support (i.e : the positions of galaxies) 
inhomogeneously distributed throughout space. Thus one has to first smooth the observed 
radial velocity field before applying no matter which inversion procedure to find the 
kinematic potential and hence the reconstructed complete velocity field. The POTENT 
inversion [2,4] has already given impressive results. However, we see two limitations 
in the smoothing scheme used in the POTENT procedure. Firstly, the size of their 
smoothing window function doesn't vary spatially, and so fictitious information is added 
in undersampled regions and a significant part of the signal is lost in oversampled regions. 
Secondly, the errors intervening during the smoothing procedure of the POTENT method 
are difficult to quantify. In this article, in section 3, we present a smoothing procedure 
based on the wavelet analysis, which permits one to optimally smooth a field sampled 
on a 3-dimensional support distributed inhomogeneously throughout space. The natural 
variation in size of the smoothing window function in our procedure permits us to realize 
a minimal or optimal smoothing procedure (i.e : without loss of information). Moreover 
we prove that our smoothed output field corresponds to a theoretical quantity (defined 
by the wavelet analysis formalism), so permitting us to quantify and control the errors 
intervening during the successive steps of our smoothing procedure. 

Finally, we analyze in a general way in section 4 the wavelet transform procedure 
involved in the reconstruction of the kinematical potential given the smoothed radial 
component of the velocity field. We show that this operation of reconstruction doesn't 
commute with the preliminary smoothing operation. This creates difficulties already at 
the 'a priori' theoretical level when attempting to compare the reconstructed velocity field 
or its associated mass density field with their correspondant quantities derived from other 
studies, since in those cases the smoothing may be effected in a different order. This 
problem is general in the sense that it is in no way unique to our smoothing procedure 
but rather depends on the global character of equation (1.1). 

2. Reconstruction of the velocity field using wavelet analysis 

2.1. The decomposition by scales 

One of the characteristic features of wavelet analysis is that it allows a simultaneous 
study of both the positional and the scaling properties of a function. Thanks to the 
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wavelet reconstruction theorem (see [6]), a square integrable function (for example the 
radial velocity field f r (x) defined at each point x of our 3-dimensional space) can be 
decomposed in a family of functions t>^(x) as follows : 

Wr (x) = (Wv r )(x) = / -^ s) (x) (2.1) 

Jo s 

The integral is performed over all scales s and f^(x) is the component of f r (x) smoothed 
on the scale s. It is in fact the spatial convolution of the radial velocity field f r (x) with 
the "reproducing kernel" K(s, x, y) (see the illustration figure la), a well-defined function 
centered on x and of spatial extension s (see [7]) namely : 

/■oo roo roo 
/ / rf 3 y K(s,x,y)v r (y) (2.2) 
-oo J -co J —CO 

v^{x) thus contains information localized "around" the scale s in frequency (scale) space. 
As the scale s decreases, a more and more detailed picture of the velocity field is thus 
available. Because of the redundant properties of the continuous wavelet transform, the 
integral involved in equation 2.1 can in practice be replaced by a discrete sum over a few 
scales (see [3,7]). Moreover, the reconstruction will necessarily be limited to a maximum 
resolution (minimum scale s c ) by the finite spatial resolution of the data distribution. 

To emphasize this latter point let us introduce two operators W Sc and W Sc acting on 
v r as follows : 

Vr = Wv r = W Sc v r + W Sc v r (2.3) 
(W. e t; r )(x) = j°° - „W(x) (2.4) 

J s c S 

(W^ r )(x) = / - „W(x) (2.5) 
Jo s 

(W Sc v r )(-x) contains the information about f r (x) at scales smaller than s c and the term 
(W Sc v r )(yL) may be considered to be a smoothed version of the radial velocity field deter- 
mined by a smoothing window function g(s c , x, y) of size s c namely : 

roo d s 

g(s c ,x,y) = / — K(s,x,y) (2.6) 

J s c S 

Note then how in practise one performs the reconstruction of v r . The v < f\'x) components 
involved in the integral over the scales in equation (2.3) are summed step by step, pro- 
ceeding from the larger cutoff scales to the smaller ones. When the desired level of spatial 
resolution is reached (for example at a given scale S), we then halt the reconstruction 
and so obtain the smoothed version of f r (x) at this scale S, (Wgi; r )(x). This procedure 
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Figure 1: For a given scale s and position x, the variations on a 2-dimensional cut of : 
(a) the "reproducing kernel" K(s,x, y), (b) the "generalized kernel" L(s, x, y). 

is particularly convenient in the case of our study for which the radial velocity field has 
to be first smoothed at an ill-defined scale s c in order to eliminate the small-scale rota- 
tional component of the cosmic velocity field. Our decomposition by scales permits us 
to add progressivly the components of the field at smaller and smaller scales, halting the 
reconstruction when the signal becomes noisy (see [7]). 

2.2. Reconstructing other quantities from the radial velocity data 

In order to derive the kinematical potential $(x) from the radial velocity field v r (x), 
we have to apply on v r the 'integral along the line-of- sight' operator P (defined in equation 
1.1). We then decompose the kinematical potential by scales as done for the velocity in 
equation (2.1), and the component & s \x.) is linked to t^(x) by the wavelet transformed 
equation (1.1): 



A practical way to reconstruct the potential is now to introduce the "generalized kernel" 
L(s, x, y) (see the illustration figure lb) by applying the operator P on the reproducing 
kernel K(s, x, y) : 




(2.7) 



L(s,x,y) 



Po^( s ,x,y) = f 1 dlK(s,lx,y) 



Jo 



(2.8) 
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It then follows from equation (2.2) that the component <3>^(x) is directly given by the 
convolution of the observed radial velocity field f r (x) by the generalized kernel L which we 
emphasize depends only on the properties of the wavelet mother function and so is known 
in principle to arbitrary accuracy, either analytically or numerically. This convolution has 
the explicit form : 



The important idea is that other quantities of physical significance may be reconstructed 
in the same way, such as the tangential components of the velocity field, the divergence 
of the field, the laplacian of the kinematical potential, etc ... . For each such quantity 
which depends linearly on the velocity field, a specific kernel function similar to L for 
the potential can be introduced that reduces the reconstruction at a given scale to a 
convolution with the radial data as in equation (2.9). For example, to the laplacian of the 
kinematical potential V 2( 3> will be associated a generalized kernel N(s, x, y) given by : 



This feature avoids in practice the amplification of the errors arising during the successive 
steps of the reconstruction procedure, especially because quantities of interest are normally 
extracted from the kinematical potential by the application of differential operators on 
$, a process that is notoriously dangerous when applied to a noisy signal. 

3. Our smoothing procedure 

3.1. The philosophy 

Peculiar radial velocities have been measured for several thousands of galaxies, and 
these have large statistical uncertainties. But the major difficulty arises from the fact that 
these galaxies with observed radial velocity are not sampled homogeneously throughout 
space. The radial velocity field v r (x) is sampled on a 3-dimensional spatial support 
(defined by the positions of galaxies) sparsely distributed, with large voids empty of in- 
formation. One thus has to first smooth the observed radial velocity data before applying 
a potential reconstruction procedure. This is done in the POTENT method [2,4] by using 
a tensor window function of large effective diameter 2500 km.s^ 1 compared to the size 
of the sample (approximatively 10000 km.s^ 1 ). One thus easily understands the impor- 
tance of the choice of the smoothing procedure and the attempts to improve upon this 
procedure (see [9]). 

The main weakness we see in the smoothing procedure used by the POTENT method 




(2.9) 




(2.10) 
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is that the effective radius of their smoothing window function is everywhere constant 
throughout the sample. Fictitious information is thus added where voids larger than the 
effective smoothing radius are present, and information is lost in oversampled regions. 
Because of the sparseness of the available data, it seems fundamental to avoid this loss 
of information. We develop below a smoothing scheme based on the use of the wavelet 
formalism which allows the radius of the smoothing window function to vary from place to 
place in the sample. In terms of wavelet vocabulary (see section 2.1), we perform the re- 
construction of the radial velocity field with a spatially variable cut-off scale s c (x), halting 
the wavelet reconstruction at a larger scale in undersampled regions and at smaller scales 
where the distribution of observed galaxies is dense. Our aim is to extract a smoothed 
velocity field containing all the information present in the catalog of the observed radial 
velocity of galaxies but no more. 

The resulting map of the cut-off scale s c (x) is wholly determined by the spatial dis- 
tribution of the 3-dimensional support (i.e : the positions {xj} i=1 jv of the N galaxies of 
the catalog). It permits us to compare our output smoothed radial velocity field derived 
from a known sparsely sampled radial velocity field {i> r ( x i)}j=i,iv, with the wavelet recon- 
struction of the true field iv(x) stopped at the cut-off scale map s c (x) : (W Sc ( x )V r )(x). 
Thus, the errors generated during the successive steps of our smoothing procedure can be 
quantified. 

3.2. The smoothed velocity field 

The observed radial velocity field we have to smooth is sampled on a 3-dimensional 
support inhomogeneously distributed throughout space. We have simulated a typical 
cosmic radial velocity field and sampled it on the support defined by the real positions 
(expressed in cartesian supergalactic coordinates) of the galaxies of the MARK II catalog 
compiled by D. Burstein (416 independent objects are present). Figure 2 shows the 
simulated radial velocity field on nine cuts passing through a cube of size 10000 km.s^ 1 
and centered on our galaxy (herein distances are expressed in km.s' 1 ; for a Hubble 
constant H of 100 km.s~ l / Mpc, the size of the sample is 100 Mpc). This non-uniform 
spatial sampling prevents us in practice from estimating integrals over space such as those 
of equation (2.2) by their associated discrete riemannian sums over the data positions. 
However, if we first restore the uniformity to the spatial support, such approximations 
can afterwards be performed with no prejudice, as long as the spatial extension s of the 
kernel involved in the convolution of equation (2.2) is greater than the elementary distance 
between neighbouring data points on the homogeneous support. Our smoothing scheme 
explores just this possibility. 

We call _E X the real 3-dimensional space wherein the spatial support {xj} i=1 7v of the 
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Figure 2: The simulated radial velocity field iv(xj) at the real positions of the MARK 
II catalog (the length of arrows is proportional to the amplitude of the field). 



Figure 3: The associated field v' r (n{x.i)) in the fictitious space. 
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N galaxies of the catalog is inhomogeneously distributed and we define by p(x) the spatial 
distribution of the support in this space. We introduce a mapping /i from this real space 
E x into a fictitious 3-dimensional space such that the image {//j = l^('X-i)}i=i,N of the 
support under the mapping ji is uniformly distributed in the space E^ : 



( — >■ E, 



det 



p(x) (3.1) 



In practice, we evaluate the mapping /i using an algorithm (this algorithm is described in 
the appendix A). 

The first step of our smoothing scheme is to associate to the set of data {tv( x i)}i=i,iv 
of the real space E x , the set {v' r (fj,i) = fr(Ai(xi))}i=i,Jv in the fictitious E^ space such that 

<(Mx)) = *v(x) (3.2) 

This operation is illustrated in figure 3 where the {f£(AO}i=i,jv are shown inside the 
normalized cube in the fictitious E^ space. 

We remark that the function v' r is now sampled on a uniform support {/Xj = /i(xi)}j =1 ^ 
in E^. It is therefore possible to define an elementary cut-off scale s^, or minimal resolution 
length, as the mean distance between 2 neighbouring points in E^ , (herein s M = 0.15) 
and to perform in the space the wavelet reconstruction W Spt v' r of v' r halted at the 
cut-off scale (see equation 2.4). The result of this operation is presented in figure 
4. Each components v^ s \fx) involved in the integral over the scales of equation 2.4 is 
estimated by replacing the spatial convolution of equation 2.2 by its associated discrete 
riemannian sum over the uniformly distributed set of points {/^}i=i,jv- Since the spatial 
extension s of the kernel K involved in the convolution is indeed always greater than 
the separation between the /Vs, this standard estimation operation can be done without 
creating fictitious information. Note that W Sfl v' r (fi) is now defined for every // of E^ and 
contains all of the information that can be extracted from the {i^(/^)}j=i,7v data in the 
fictitious E^ space. 

The last step of our smoothing procedure is to return to the real space E x through 
the inverse mapping Our output smoothed velocity field (Mf r )(x) is finally the field 
corresponding to W Sti v' r in the real space E x namely : 

(M Wr )(x) = (W^XXMx)) (3-3) 

This operation is illustrated figure 5. Since the mapping /i establishes a one-to-one cor- 
respondance between points of the real space E x and those of the fictitious E^, our 
output smoothed radial velocity field (Mv r )(x) contains the maximum of information 
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Figure 4: The wavelet reconstruction (W Sll v' r )(/j,) of v' r halted at the cut-off scale s M . 

which can be extracted from the {fr (/■**) }i=i,w and thus from the observed radial veloci- 
ties {v r (xi)} i= i t N of the N galaxies in the catalogue. In this way our smoothing procedure 
is minimal (no loss of information). 

3.3. Correspondence with a theoretical quantity 

Thanks to the scale resolution properties of the wavelet transform, our output smoothed 
radial velocity field (Mv r )(x) can be expressed as a theoretical quantity. We prove in ap- 
pendix B that, so long as the mapping ji satisfies a validity condition (see below equation 
3.5), the following equality holds : 



Thus our output smoothed radial velocity field may be identified with the wavelet recon- 
struction of v r halted at the inverse map of the cut-off scale s c (x) which varies with x. 
We show in figure 6 the cut-off scale inverse map corresponding to the spatial distribution 
previously presented in figure 2. The value of s c (x) is derived from the jacobian associated 
with the mapping \x (see appendix A). The lower the density at the position x, the larger 
is the corresponding cut-off scale. 

We exhibit in figure 7 the wavelet reconstruction of the previously simulated radial 
velocity field halted at the cut-off scale s c (x). We notice that even if the main features 





s 



(3.4) 
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Figure 5: Our output smoothed radial velocity field (Mv r )(x) (in the real space E x ). 



remain, our smoothed radial velocity field (Mv r )(x) of figure 5 differs in detail from 
(W Sc ( x )V r )(x). The reason for this discrepancy is that the mapping // doesn't in fact sat- 
isfy the validity condition (appendix B) which stipulates that for every x and vector h 



if ||h|| < s c (x) , 



dx k 



.[h] 



det 



dxk 



1/3 



x llhl 



(3.5) 



or in other words that the mapping fi is locally equivalent to a rotation-dilation transfor- 
mation (see appendix B). 

However, this discrepancy between our output smoothed velocity field (M« r )(x) and 
the wavelet reconstruction (W / Sc ( x )f r )(x) of the true radial velocity field f r (x) halted at 
the cut-off scale map s c (x) can be quantified by analysing the spatial distribution of the 
N galaxies of the sample. Thus at each position x in the sampled volume, we can define 
a "relative sampling error" E(x) computed as follows : 



E 2 ( 



Ei=i [^(sA»M x )>My»)) - 1 /p(yi) K ( s c{*),x,yi)Y 



(3.6) 



This 'relative sampling error map' E(x) gives at each position x the relative error between 
(Mtv)(x) and (W / Sc ( x )f r )(x) as a fraction of the amplitude of the true radial velocity field 
IV (x). The relative sampling error map -E'(x) for the spatial distribution of the galaxies 
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Figure 6: Isocontours of the cut-off scale map s c (x) (Normal, heavy and dotted contours 
are resp., 750, 1500 and 2250 km.s^ 1 ). 



Figure 7: The wavelet reconstruction (W Sc ( x )iv)(x) of the simulated radial velocity field 
IV (x) halted at the cut-off scale map s c (x). 
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Figure 8: Isocontours of the relative sampling errors map E(x) (Normal, heavy and 
dotted contours are resp., 15%, 25% and 50%). 

of the MARK II catalogue is shown in figure 8. In practice we can discard the regions of 
the sample where the sampling errors are too high. 

4. The kinematical potential $(x) = (Pt> r )(x) 

We have shown in section 3.3 that, in the regions of the sampled volume where the 
validity condition is verified (i.e : the domains where the relative sampling errors -E'(x) 
are low), our smoothed radial velocity field (Mw r )(x) can be identified with (W Sc ^v r )(x). 
As we suggested in the introduction, there are good reasons to believe that the cosmic 
velocity field is irrotational up to an ill-defined scale and thus that a kinematic potential 
$(x) can be inferred by integrating the cosmic radial velocity field tv(x) along the line- 
of-sight (equation 1.1). It turns out that the ill-defined scale of irrotationality is generally 
presumed to be smaller than the cut-off scales s c (x) which limit the spatial resolution of 
the accessible cosmic radial velocity field (W Sc ( x )tv) (x). It becomes thus possible to extract 
the velocity potential $(x) by applying the 'integral-along-the-line-of-sight' operator P 
either to (W / Sc ( x )f r )(x) or to the smoothed radial velocity field (Mt> r )(x). 

Unfortunately, this operator P has a non-local character. Hence, the potential derived 
from the smoothed radial velocity field PoW Sc ( x )V r differs from the smoothed potential of 
the velocity field W / Sc ( x ) oPv r (and of course from the non-smoothed kinematical potential 
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Figure 9: Isocontours of the potential of the smoothed radial velocity field (P o 
W / s c (x)^r)(x) (the contour spacing is 2.5 10 5 (km.s' 1 ) 2 , negative contours are dotted, heavy 
contour is 0). 



Figure 10: Isocontours of the smoothed simulated potential (W Sc ^oPv r )(yi) (the contour 
spacing is 2.5 10 5 (/cm.s -1 ) 2 , negative contours are dotted, heavy contour is 0). 
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$ = Pv r ) (see appendix C). We illustrate this discrepancy by plotting in figure 9 the 
potential derived from (W Sc ( x )i>r)(x) and in figure 10 the smoothed simulated potential 
(W / Sc ( x )$)(x). We want to emphasize that this behaviour is not due to the way we smooth 
the radial velocity field but is intrinsically linked to the nature of the operator P. 

The point has its importance since the kinematical potential (or the reconstructed 
3-dimensional velocity field) derived from catalogues of radial peculiar velocities is often 
considered as reference data for other studies (for example Saunders et al. (1991), Dekel 
et al. (1993)). As an example note that, above the scale of irrotationality, the total 
mass density perturbation field 5 t (x) is linked to the kinematical potential <3>(x) through 
the Poisson equation (S t (x.) oc V 2 $(x)). Adopting the hypothesis that the luminous 
matter traces the mass (5 t (x) oc 5/(x)), a smoothed total mass density field may be 
extracted from observations of the spatial distribution of galaxies. Unfortunately for the 
intercomparison of such independent results : 

(W Sc 5 t )(x) = {V\W Sc o P Wr ))(x) + (V 2 (P o W Sc v r ))(x) (4.1) 

One thus should be very cautious when comparing the density or kinematical potential 
derived from observed radial peculiar velocity catalogues with similar quantities obtained 
from other studies, such as those based on number galaxy counts. 

5. Conclusions 

We have presented a method, based on the properties of the wavelet transforms, for 
smoothing a field sampled on a support inhomogeneously distributed throughout the space 
of interest. Our smoothing scheme is minimal (no loss of information) and our output 
smoothed field can be identified with a well-defined theoretical quantity, as long as the 
spatial support of the field satisfies certain criteria. We expect this technique to be of 
quite general importance. The particular application of this smoothing scheme to the 
observed cosmic radial velocity field discussed above, reveals some apparently previously 
unrecognized limitations concerning the reconstruction of the kinematic potential from 
the smoothed radial velocity field. Indeed, we prove that this potential will generally 
differ a priori from a similar quantity obtained from other cosmological studies. The 
existence of this kind of error now has to be taken into account. 
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In practice, we evaluate the mapping // using the following algorithm. We first choose 
a direction (the Xi-axis for example). Along this direction, we sort the N points of the 
sample by the increasing value of their xi-coordinate. To the object of rank i, we then 
affect the value ^ = i/N to its /^-coordinate, the first component of the mapping /i. 
This procedure ensures us that the set of the /ii-coordinates of the sampled objects are 
uniformly distributed along the /ii-axis in the fictitious space E^. We now divide the 
cube in slices perpendicular to the rri-axis (7 rri-slices for this sample) such that each slice 
contains approximatively the same number N\ of data points. The non-uniform sampling 
in the E x space implies that each slice have its own width. Nevertheless, the widths of 
these slices become equal along the /xi-axis in the space E^. 

The next step of our algorithm is to repeat this construction for each Xi-slice, choosing 
the x 2 -axis as the sorting direction. Inside each xi-slice, the object of rank j has a 
A*2 = j/Ni. We thus obtain /^-coordinates uniformly distributed inside each a^-slice, 
and so for the set of all slices in the fictitious E^ space. The Xi-slices are afterwards 
divided again in slices perpendicular to the X2-axis (7 for this sample) containing the 
same number N 2 of objects and each xiX2-slice is sorted along the X3-axis (for an object 
of rank k, = k/N^). N 2 is 8 or 9 in our case, depending on the xiX2-slice treated. 

At the output of this algorithm, we then obtain a set of N points {/ij = l^('X-i)}i=i,N 
uniformly distributed inside a normalized cube in the fictitious space E^. A continuous 
version of the /i mapping as well as its inverse mapping /i -1 can be computed, if necessary, 
by using a linear interpolation on the fi\, n 2 and fi^ components of the transformation /i. 

In order to evaluate quantities such as the cut-off scale map s c (x) (section 3.3.), the 
spatial density distribution of the support p(x) in the real space _E X has to be evaluated. 
We then first estimate at the position x the 9 terms (dfij/dxk)(x-) (j and k vary from 1 
to 3) of the transformation coordinate matrix associated with the mapping /i. We then 
derive from this matrix its jacobian J M (x) and finally adopt for the density distribution 
p(x) a smoothed version of J M (x). 

Appendix B : Comparison between (Mv r )(x) and (W / Sc ( x )f r )(x) 

Our goal here is to identify our output smoothed radial velocity field (Miv)(x) with 
a theoretical quantity derived from the true radial velocity field t> r (x). Our smoothing 
procedure ensures us that the following equality holds (see equation 3.3) : 

(Mtv)(x) = (W.X)(/i(x)) (B.l) 

In the fictitious E^ space, (W Sll v' r .)([i) is the wavelet reconstruction of the field v r (fx) halted 
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at the cut-off scale s M . It satisfies (see equations 2.3, 2.4 and 2.5) : 



(B.2) 



Now, Let us focus on the measurable component (W Sfl i^.)(/z(x)). In the fictitious 
space, it can be expressed as follows (see equations 2.2 and 2.5) : 



l-Sa A a f 

(WW r )(^)) = - d 3 rv' r (r)K(s,^),r) 
Jo s Je u 



(B.3) 



Recalling equation 3.2 and introducing y such that r = /i(y), this equation becomes in 
the real space _E X : 



(w^;)(Mx)) = 



ds 
s 



P(y) d 3 yv r (y) K(s, /i(x), //(y)) 



(BA) 



At this stage, we have to use the specific properties of the reproducing kernel K. This 
kernel K can be derived from an "analysing" function k of unit spatial extension due to 
the following symetry (see [7]) : 



{B.5) 



Hence, for every real number a, the kernel satisfies : 

1 



K(s, ax, ay) 



a u 



(B.6) 



Keeping this in mind, we develop the mapping fi to the first order of its vectorial taylor 
series with respect to x or y. The difference /i(x) — /i(y) involved in the kernel K 
intervening in equation B.4 can then be approximated : 



IIM x ) - My) 



dx k 



(x) 



[x-yj 



djij 
dx k 



y) 



[x-yj 



(B.7) 



This approximation is in practice performed inside regions of the space not larger than 
the cut-off scale (these regions are in fact defined by the spatial extension of the kernel 
K involved in the spatial convolution of the equation B.4). 

The next step is to make the crucial assumption that the mapping \i can be locally iden- 
tified with a rotation-dilation transformation, or in other words that for ||h|| < s^j p 1 ^ 3 (x) 



d/ij 
dx k 



x) 



det 



dx k 



1/3 



x ||h|| = p(x) 1/3 x ||h| 



[Bl 
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This equation is certainly not valid for every point x of the sampled volume. Nevertheless, 
in the regions of space where this validity condition does hold, the following approximation 
is valid (see equations B.5, B.6, B.7 and B.8) : 

if s < Sfl , K( S ,/i(x),/i(y)) » -l_X( s /p( x )V3 )X) y) (5.9) 

Then we can define a cut-off scale map s c (x) in the E x space : 

and it turns out that if we introduce S such that S = s/p(x.) 1 ^ 3 , equations B.4, B.9 and 
B.10 give : 

(WW r )(^)) » / 4 / d 3 yv r (y)K(S,x,y) = (W°™v r )(x) (B.ll) 
Jo b Je x 

Finally, because of equations 2.3, B.l, B.2 and B.ll, the following equality holds in the 
regions of the sampled volume where the validity condition of equation B.8 is satisfied : 

(Mu r )(x) = (W Sc{x) v r )(x) (5.12) 



Thus, our output smoothed radial velocity field (Mt> r )(x) can be compared to the 
wavelet reconstruction (W Sc ( x )i; r )(x) of the real radial velocity field f r (x) halted at the 
cut-off scale map s c (x). In practice we evaluate the "degree of approximation" of the 
validity condition by comparing the 2 sides of the equation B.9 integrated over the space 
at the fixed scale s M . This quantity i?(x) is defined as follows : 

E 2 ( x E£i [K(s„,fi(^),n(yj)) - l/p(yQX(g c (x),x, yi )] 2 

This relative sampling errors map E(pc) gives at each position x the relative error between 
(M« r )(x) and (W Sc ( x )i; r )(x) with respect to the amplitude of the true radial velocity field 
u r (x). 

Appendix C : P o W Sc{x) v r ^ W Sc(x) o Pv r 

We prove in this appendix that the 'integral-along-the-line-of-sight operator' P doesn't 
commute with the necessary smoothing operation that must first be performed on the 
radial velocity field iv(x). This behaviour doesn't depend on the particular nature of 
our preliminary smoothing, rather it remains true more generally even if for example the 
effective radius of the smoothing window function is constant throughout the sampled 
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volume as in POTENT. 

Let us express the potential derived from the smoothed radial velocity field P o W Sc v r 
in terms of the smoothing window function introduced equation 2.6 (herein , s c is chosen 
constant throughout the space since otherwise our argument is even stronger) : 

(PoW Sc v r )(x) = ^ dl f d 3 yv r (y)g(s c ,lx,y) (C.l) 
Jo Je x 

(PoW Sc D r ) (x) is the integral along the line-of-sight of the field v r smoothed with a window 
function of constant effective radius s c from O to x. 

In the same way, let us express the smoothed version at the scale s c of the potential 
(W Sc oPv r )(x) : 

(W Sc oPv r )(x) = I d 3 y ^ dlv r (ly) g(s c ,x,y) (C.2) 

JE X JO 

Because of the specific properties of the wavelet smoothing window function, it turns out 
that (see equation 2.6, B.5 and B.6) : 

g(s c ,x,y) = I 3 g(ls c ,lx,ly) (C.3) 

Hence equation C.2 can be re-expressed, if we introduce Y such that Y — ly : 

(W Sc oPv r )(-K) = f 1 dl [ d 3 Yv r (Y)g(ls c ,lx,Y) (C.4) 
Jo Je x 

(W Sc o Pt> r )(x) is then the integral along the line-of-sight of the field v r smoothed with a 
window function with a variable effective radius ls c all along the line-of-sight, from at 
the observator position O to s c in x. It thus follows that : 

PoW Sc v r + W Sc oPv r (C.5) 

as was required. 
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